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Abstract. As has been recently realized, a certain two-point function A^v - and its associated form factors G and L - play a 
prominent role in the PT-BFM formulation of the Schwinger-Dyson equations used to study gauge-invariantly the gluon and 
ghost propagators. After showing that in the (background) Landau gauge A^y fully constrains the QCD ghost sector, we show 
that G coincides with the Kugo-Ojima function u, whose infrared behavior has traditionally served as the standard criterion for 
the realization of the Kugo-Ojima confinement mechanism. The determination of the behavior of G for all momenta through 
a combination of the available lattice data on the gluon and ghost propagators, as well as the dynamical equation G satisfies, 
will be then discussed. In particular we will show that in the deep infrared the function deviates considerably from the value 
associated with the realization of the Kugo-Ojima confinement scenario; the dependence on the renormalization point of u, 
and especially of its value at = 0, will be also briefly discussed. 
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Over the last few years large volumes ab-initio lat- 
tice gauge theory computations have established beyond 
any reasonable doubt that the gluon propagator and the 
ghost dressing function of pure Yang-Mills theories in 
the Landau gauge saturates in the deep infrared (IR) 
at a finite, non-vanishing value both for SU(2) [1] and 
SU(3) [2] gauge groups. Specifically choosing the 
Landau gauge and defining the gluon propagator cofactor 
A, and the ghost dressing function F as 



(1) 



where Pu y (?) = gfiv ~ q^iQv /q^ is the transverse projec- 
tor, and D{q^) the ghost propagator, the aforementioned 
lattice results tell us that (in Euclidean space) 



A-'(0)>0, and F(0) > 0. 



(2) 



In the continuum formulation, the only way of ob- 
taining these so-called massive solutions in a gauge in- 
variant way and without breaking (either explicitly or 
softly) the BRST symmetry of the original Yang-Mills 
action, is within the PT-BFM framework [3], where a 
truncation scheme that respects gauge invariance at ev- 
ery level of the dressed-loop expansion has been devel- 
oped in [4] using the systematic rearrangement of the en- 
tire Schwinger-Dyson series allowed by the pinch tech- 
nique [5, 6, 7]. 

In the PT-BFM construction one studies the PT-BFM 
propagator A which is related to the conventional prop- 
agator of Eq. (1) through the background-quantum iden- 

tity [8] 

A{q^) = [l+G{qYHq^), (3) 
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FIGURE 1. The auxiliary functions Apy and Hav appearing 
in the PT-BFM framework. 



where G is the g^y form factor appearing in the Lorentz 
decomposition of the auxiliary Green's function A^y 
defined as (see Fig. 1) 

^tiv{q) = s'-Ca [D{k + q)A"^{k)Hav{k,q) 

J k 



g,Mq') + '-^L{q^ 



(4) 



The function //^y appearing above (Fig. 1 again) is in 
fact a familiar object, since it appears in the all-order 
Slavnov-Taylor identity satisfied by the standard three- 
gluon vertex. It is also related to the full gluon-ghost 
vertex by the identity q^Hnv{k,q) = —ir^{k,q); at 

tree-level, i/^y' = ig^v and r''^\k,q) r^{k,q) = -q^. 

The G and L form factors play a prominent role in the 
(background) Landau gauge, where the presence of an 
extra local functional equation (the so-called antighost 
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FIGURE 2. Left panel: The form factor —G{q^) determined from Eq. (24) at different renormalization points /i through the 
procedure described in the text. Right panel: Same as in the previous panel but this time for the L{q^) form factor. 



equation) implies the identity [9] 



F-\q^) = \+G{q^)+L{q^) 



(5) 



a relation that is valid also in the conventional Landau 
gauge [10]. Since, under very general conditions on the 
gluon and ghost propagators, L{q-^) — > when — > 
one has the IR relation F^'(O) = 1 +G(0). Thus, we 
see that a divergent - or enhanced [1 1] - dressing func- 
tion requires the condition G(0) = — 1 . To the practition- 
ers, this latter condition will look suspiciously similar to 
the Kugo-Ojima (KO) confinement criterion which de- 
mands (as a necessary condition for confinement through 
the so-called quartet mechanism) that a certain func- 
tion u[q^) (the KO function) acquires the IR value 
m(0) = —1 [13]. Indeed, it is possible to show that G is 
nothing but the KO function [9, 12] 



J-\ — 



)^G{q'). 



(6) 



Therefore the form factor G encodes practically all rele- 
vant information on the IR dynamics of the ghost sector, 
and, at least partially, the gluon sector as well [through 
the identity (3)]. 

Approximating the three point functions and //uv 
with their tree-level value (the first approximation on the 
ghost-gluon vertex being supported by lattice studies), 
the dynamical equations satisfied by G and L read 



G{q^) = 



L{q^) = 



(k-qf 



1-3 



(k-q) 



k^q^ 



A{k)D{k + q), 
21 

Aik)D{k + q).i7) 



Then, since within our approximation scheme, in the 
equations above only A, F (through the ghost propagator 
D) and g appears, in order to determine the behavior of G 



and L one can fully exploit the available lattice data on A 
and F , through the following general strategy [12, 14]. 
One starts by using the lattice gluon propagator as an in- 
put for the ghost SDE; then solves for the ghost dress- 
ing function, tuning the coupling constant g such that 
the solution gives the best possible approximation to the 
lattice results. Obviously one must check that the cou- 
pling so obtained (at the renormalization scale used for 
the computation) is fully consistent with known perturba- 
tive results (obtained in the MOM scheme, which is the 
scheme used in our computations); this is indeed what 
happens [12, 14]. At this point the one has the three 
building blocks A, F and g fully determined, and can start 
analyzing other quantities constructed from them such as 
the G and L form factors above [12] or the renormaliza- 
tion group invariant effective charge [14, 15]. 

Before solving numerically the equations (7), there is 
one last issue that needs to be addressed. Specifically, 
one needs to identify a renormalization procedure for G 
and L that does not break the identity (5), which, due to 
its BRST origin, should not be deformed (within the PT- 
BFM scheme) by the renormalization process. Note in 
fact that Eq. (5) constrains the cutoff-dependence of the 
unrenormalized quantities involved; specifically, denot- 
ing by Zc the ghost wave-function renormalization con- 
stant (Z^Fq"' = F^^) and with Za the (yet unspecified) 
renormalization constant of the function A^v(^), with 

Za^'"' + Aq"] ^ + A^", one finds that (5) is pre- 
served iff Za =Zc [15]; as a result, one finds the relation 

Z,(A^M^)[l+Go(<?^A2)+Lo(^^A2)] 

= l+Giq\^^)+L{q\^^). (8) 

Imposing then the renormalization condition F{pL^) = 1, 
going to Euclidean space, setting q^ =x,k^ =y and a, = 
g^ /An, and implementing the standard angular approxi- 



mation, one finds the renormalized equations [12, 15] 



1+G(x) = Z,- 



16;r 



F{x) 



jyy(3 + y]Aiy)Fiy) 



L{x) = 



12% 



F{x) 



dy 



x^ JO 
A{y)F{y) 



dyy^Mj) 



y 



(9) 



Notice that L is finite, as expected from power counting; 
in addition, we see (by means of the change of variables 
y = zx) that if A and F are IR finite, then L(0) = 0, as 
mentioned before (notice however that the same result is 
obtained for scaling solutions [11], where A(3;) ~ y" and 
F [x) ^ x^, provided that a + b> — 1 ). 

At this point, all necessary ingredients for determining 
the functions G and L are available. Substituting them 
into the corresponding equations given in (9), we ob- 
tain the solutions shown in Fig. 2, where we see that L 
is subdominant, and that indeed it vanishes in the deep 
IR. Also, the jx dependence of the KO function and the 
KO parameter are clearly shown; for the range of renor- 
malization points ^ chosen, the KO function saturates in 
the deep IR to the value G(0) k, —0.6 which deviates ir- 
remediably from the value G(0) = — 1 required for the 
realization of the KO confinement scenario. 

The curves plotted for G should finally be compared 
with those obtained on the lattice in [16], where the KO 
function u was studied in terms of Monte Carlo averages, 
and its asymptotic behavior inferred from the identity 
(5). Even though in [16] the extrapolation towards the 
zero limit was problematic, due to a lack of knowledge 
of the function L (our analysis does not suffer from such 
a limitation, given that L is completely determined by its 
own equation) and that, for essentially the same reason, 
the renormalization procedures employed are different, 
we clearly see in Fig. 3 the same behavior emerging, and 
in particular the saturation of the KO function in the deep 
IR to a value sensibly different from the critical — 1 . 

In conclusion, we have shown that the massive gluon 
propagator A and ghost dressing function F found as 
solutions to the SDE and confirmed by all large vol- 
ume lattice simulations up to now, do not support a 
confinement scenario based on the original KO mecha- 
nism/criterion. However, due to the equality (6) between 
the KO function u and the auxiliary function G, and the 
central role that the latter function has in the PT-BFM 
scheme, it would still be very interesting to carry a 
thorough study of such function on large volume lattices 
(for different space-time dimensions and gauge groups). 
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FIGURE 3. Comparison between our results and direct lat- 
tice calculation of the KO function —uicp-) at /i = 4 GeV. 
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